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ABSTRACT 

yA/ms. We report the first observation of multiple-periodic propagating disturbances along a fan-like coronal structure simultaneously 
detected in both intensity and Doppler shift in the Fexii 195 A line with the EUV Imaging Spectrometer (EIS) onboard Hinode. A 
new application of coronal seismology is provided based on this observation. 

Methods. We analyzed the EIS sit-and-stare mode observation of oscillations using the running diff'erence and wavelet techniques. 
Results. Two harmonics with periods of 12 and 25 min are detected. We measured the Doppler shift amplitude of 1-2 km s~\ the 
relative intensity amplitude of 3%-5% and the apparent propagation speed of 100-120 km s~^ 

Conclusions. The amplitude relationship between intensity and Doppler shift oscillations provides convincing evidence that these 
propagating features are a manifestation of slow magnetoacoustic waves. Detection lengths (over which the waves are visible) of the 
25 min wave are about 70-90 Mm, much longer than those of the 5 min wave previously detected by TRACE. This diff'erence may be 
explained by the dependence of damping length on the wave period for thermal conduction. Based on a linear wave theory, we derive 
an inclination of the magnetic field to the line-of-sight about 59+8°, a true propagation speed of 128+25 km s~^ and a temperature of 
0.7+0.3 MK near the loop's footpoint from our measurements. 
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1. Introduction 

Quasi-periodic propagating intensity disturbances in fanlike 
large co ronal loops were first obse rved with SOHO/EIT 195 
A data (iBerghmans & Clette1ll999|), and then we re confirmed 
with TRACE 171 A data (e.g. iDe Moortel et al. II2OO Q). These 
disturbances with a propagation speed on the order of 150 
km s~^ have been interpreted as slow magnetoacoustic waves 
' (iNakariakov et al. II2OQQ) . Detection of these waves is important 
, for our understanding of the energy balance in the outer solar at- 
mosphere and also hi ghly valu able for coronal seismology (e.g. 
reviews by Nakariako v & Verw ichte 2005). 

Many authors have found evidence of 3 min and 5 min 
oscillations propagating through the c hromosphere and transi- 
tion region into the lower corona je^g. lDe Moortel etal.ll2002t 
iMarsh et al. ll2003l: IWang et al. "2009"). supporting their origin in 
the p hotospheric /?-modes by wave leakage (De Pontieu et al. 
I2005h . In addition, there were also a few reports of similar 
propagating waves with long periods (10-15 min) observed 
in coronal loops by EIT (IBerghmans & Clette I [T999h. TRACE 
("Mcintosh et^ l2QQ8h . and STEREO/EUVI (Marsh etaD 
i2009), but their origin is still not clear. Long-period (10-30 
min) oscillations wer e also observed in polar coronal holes (e.g. 
iBaneriee et ani2009h . 

In this study, we detect two long-period (12 and 25 min) 
harmonics in propagating slow magnetoacoustic waves in fan- 
like coronal loops with observations from the EUV Imaging 
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Spectrometer (EIS) onboard Hinode. The EIS instrument 
with high sensitivity and spectral resolution (ICulhane et al. I 
I2QQ7D can achieve relative Doppler- shi ft measurements with 
an accuracy of less than 1 km s~^ (iMariska et al. I I2QQ8I: 
Ivan Doorsselaere et al. I l2QQ8l: lErdelvi & TarovanI l2008b . pro- 
viding us a good opportunity for detecting oscillations and waves 
in solar atmosphere. 



2. Observation 

EIS has both imaging and spectroscopic capabilities. Its spectro- 
scopic mode can operate in a rastering mode or a sit-and-stare 
mode. In the latter case, the V or slit is placed at a fixed lo- 
cation on the Sun, and repeated exposures are obtained while the 
Hinode spacecraft tracks the solar rotation. 

The observations were obtained on 2007 February 1 in AR 
10940. An EIS spectroheliogram was taken from 00:12 to 01:22 
UT with the V slit and covering a 256'' x 256'' region. A sit- 
and-stare observation within the region was taken from 01:32 to 
07:32 UT with the 1" x 512" slit and an exposure time of 60 s. 
The data cover 17 spectral windows with a cadence of about 
62 s. This paper presents the result mainly for the Fexii A 195 
line. The raw data were first processed using the standard rou- 
tine eis.prep in Solar Software (SSW). The Fexii 195.12 A line 
was then fitted with a single Gaussian profile, providing the to- 
tal intensity, Doppler shift, and line width. The eflTect of a weak 
blended line, 195.18 A (lYoung et al. ll2Q09h . is negligible since 
only relative Doppler shift variations are concerned in our study. 
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(b) EIS Fe XII 00:12 - 01 :22 UT (c) Doppler Velocity 
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Fig. 1. (a) The TRACE A195 bandpass image, (b) The intensity map in the Fe xii /1 195. 12 line from EIS. The raster observation was 
taken from 00:12 to 01:22 UT on 2007 February 1. (c) The Doppler velocity measurements. The red color represents the redshift 
and the blue color the blueshift with a scale range from -20 km s"^ to +20 km s"^ . The vertical line in each plot shows the position 
of the EIS 1 slit for the sit- and- stare observation. The short horizontal lines on the slit mark the positions where the oscillations are 
analyzed. 



Hinode is known to have instrumental jitter in both the 
X and y directions. We estimated drifts of the pointing us- 
ing t he SSW routine eis Jitter from the satellite alignment data 
(IShimizu et al. Il2007i) , and find that the displacements of the EIS 
pointing are within 3'' during the observation. For the sit-and- 
stare data analyzed below, the pointing drifts in the }^-direction 
have been corrected. The drifts in the x-direction have little ef- 
fect on our result since they are slow and orbit-related with an 
amplitude less than and a period of about 90 min. 

3. Results 

A part of TRACE 195 A image recorded at 03:15 UT shows 
a typical fan-like set of diverging coronal loops, located at the 
northern edge of NOAA active region 10940 (Fig.[It). The EIS 
spectroheliogram taken at 00:12 UT only covers part of this 
structure (Fig.[T]3). The slit used for the sit-and-stare observation 
is almost aligned along a loop. The Doppler shift map in EIS 
Fexii /1 195 shows that the eastern part of footpoints of the fan- 
like loops is dominated by redshifts, while the western part of 
footpoints is dominated by blueshifts (Fig.[T]:). The sit-and-stare 
mode observations show weak blueshifts less than 7 km s"^ at a 
region of j=90"-120". 

Time-distance maps in intensity and Doppler shift for the 
EIS Fe XII line are shown in Fig.[2l To enhance weak propagating 
intensity disturbances, the time series of relative intensity were 
normalized to the background trend at each position along the 
slit by [I(t, yi) - hg(t, ydVhgit, yd (Fig.[2^), where hg(t, yd is the 
background trend aty = y^ taken as a 20 min average smoothing 
of I(t,yi), while the 30 min average smoothing gives almost no 
diff'erence. The background trend was also subtracted from the 
time series of Doppler shift with a similar method (Fig.[2ll). The 
upwardly propagating disturbances can be seen in both intensity 
and Doppler shift as weak bright and dark bands with the posi- 
tive gradient. The propagating disturbances are slightly aff'ected 
by a small bright point at j= 130 ' -140 '. The footpoint of the 
loop (at j=80''-100'') appears to be characterized by the rela- 
tively high-frequency features, while the upper part is dominated 
by the relatively low-frequency features. In other words, the 



lower frequency component appears to propagate over a longer 
distance than the higher frequency one. These time- variable fea- 
tures can be enhanced in a running diff'erence image, which was 
created by subtracting from each frame a frame taken some ex- 
posures later. The middle and right panels of Fig.[2l show the 
5 min and 10 min running diff'erence images, which clearly re- 
veal the presence of multiple-period modes in the propagating 
disturbances. The propagation velocity transverse to the line of 
sight is estimated as the gradient of the lines outlining the evi- 
dent slanting bands in the running diff'erence images (see two ex- 
amples shown in Figs.[2j3 and[2t). For nine propagation features 
selected for intensity and Doppler shift separately, we measured 
the projected propagation speed of 118+34 km s"^ for intensity 
and 105+25 km s"^ for Doppler shift. In addition, we also exam- 
ined the coordinated TRACE data, but no corresponding propa- 
gating features are found in the constructed time-distance image 
along the EIS slit (see Fig. lA.ll online). This may come from the 
lower sensitivity of TRACE compared to EIS. 

The periods of the distu rbances were invest i gated using 
a wavelet analysis method dTorrence & Compol Il998h . The 
Morlet wavelet was chosen for the convolution of time series. 
The intensity and Doppler shift oscillations at three positions 
(y=l50 " , no ', and 60 ' -65 ', denoted as oscillations 1, 2, and 3) 
were selected for analysis. Oscillation 1 is located at the position 
of a loop with a projected distance of about 50 Mm away from 
the footpoint, near which oscillation 2 is measured. Oscillation 
3 was selected for comparison, located at another loop intersect- 
ing the slit (Fig. [It). The time profiles used for the analysis of 
oscillation 3 were obtained by averaging over six pixels along 
the slit because the oscillation region varied slightly in the y- 
direction during the observation. In practice, we first subtracted 
the background trend for both intensity and Doppler shift, which 
was obtained by temporally smoothing the series. Since oscilla- 
tion 1 is dominated by the lower- frequency wave component, the 
background is taken as a 20 min average smoothing. To empha- 
size the shorter- frequency wave component in oscillations 2 and 
3, the background is chosen as a 10 min average smoothing. 

The evolution of detrended relative intensity and Doppler 
shift for three oscillations is shown in panels (a) of Figs. [3] and 



Wang et aL: Propagating low-frequency slow magnetoacoustic waves in coronal loops 



3 



Fe XII Intensity I (5 min r ii [ ii J i ^ ^ i - i - I (10 mm running difference) 




Fig. 2. The upwardly-propagating waves in coronal loops observed by Hinode/EIS. (a) Time series of relative intensity along the 
slit in the Fexii /1195.12 line, (b) and (c) The 5 min and 10 min running difference maps for intensity, (d) Time series of Doppler 
shift. Here the white color indicates the blueshift and the black color indicates the redshift. (e) and (f): The same as (b) and (c) but 
for Doppler shift. The horizontal lines in (a) and (d) mark three positions (y=l50'\ llO ', and 60 ' -65' ), where the time series of 
oscillations are analyzed (see Figs.[3]and|4l and Fig. lA.2l (online only). The lines in (b) and (e) outline the propagation features along 
the slit to demonstrate measurements of the propagation speed. 



m and of Fig. lA.2l (online only). There is a good in-phase re- 
lationship between the intensity and Doppler velocity, which is 
consistent with the signature as predicted by linear MHD wave 
theory for an upwardly propagating slow magnetoacoustic wave 
in coronal loops. The wavelet spectrum and the global wavelet 
spectrum are shown in panels (b)-(e) of Figs.[3] and (H and of 
Fig. lA.2l (online only). Cross-hatched regions in the wavelet 
spectrum indicate the "cone of influence", where edge efl'ects be- 
come important due to the finite length of time series. The global 
wavelet spectrum is the average of the wavelet power over time 
at each oscillation period. 

For oscillation 1 the wavelet spectra for both intensity and 
Doppler shift shows strong power at the period in a range of 
20-30 min over a duration of 6 cycles. From the global wavelet 
spectra, the oscillation period for Doppler shift is measured to 
be 25.6+6.3 min and the period for intensity is 26.2+4.9 min. 
The uncertainties are taken as the half FWHM of the peak. The 
average maximum amplitude for Doppler shift is 1 .6+0.5 km s"^ 
and for relative intensity is (4.9+1.5)%, where the errors are the 
standard deviation of absolute peak values. As a comparison, we 
also examined the case with a subtraction of the 30 min smooth- 
ing background (see Fig. |A3] online), and find that the power of 
the 26 min band shows almost no change, while the power of the 
88 min band for Doppler shift rises above the 99% confidence 
level, which is caused by the orbital eff"ect. 

Oscillation 2 clearly shows that most of the wavelet power 
is concentrated within two period bands ranging from 8-16 min 
and from 20-30 min, indicating that the waves near the footpoint 
consist of multiple harmonic modes. The two periods are mea- 
sured to be 11.9+2.8 min and 25.1+7.0 min for Doppler shift, 
and 13.5+4.7 min and 26.2+5.0 min for intensity. The average 
maximum amplitude for Doppler shift is 1.3+0.3 km s"^ and for 
intensity is (3.8+0.6)%. Oscillation 3 shows most of the wavelet 
powers for Doppler shift and intensity concentrated in a range 
of 8-14 min. The measurements of periods and amplitudes for 
oscillations 1-3 are summarized in Table l A. ll (online only). 

We show a new application of coronal seismology based on 
interpreting our observation in terms of slow magnetoacoustic 



waves. Assuming that the coronal field, along which the waves 
propagate, has an inclination, 0, to the line of sight, we have 



Vn 



sm0 = 



(1) 



where Vp is the projected propagation speed of the waves and 
Cs the sound speed, which is a good approximation to the phase 
speed of slow magnetoacoustic waves under a typical coronal 
condition of low-/? plasma. The linearized continuity equation 
gives pVPo = y /Cs, where po is the undisturbed loop density, 
and p' and V the density and velocity perturbations, respec- 
tively. Considering r /Iq ~ 2p7po and cos0 = V\\/V\ where 
r/Io is the relative amplitude for intensity disturbances and V\\ 
the Doppler shift amplitude, we obtain 



2Vu 



COS0 = 



(r/k)cs 

From Eqs. ([U and (|2]) we derive the inclination angle. 



tan 



1 

2^^ 



^1 



(2) 



(3) 



Note that Eq. 0]) also holds in the condition with a steady flow in 
the loop. We take 1 10+30 km s"^ , the mean of the projected 
propagation speeds measured for intensity and Doppler shift. By 
assuming that the inclination of the coronal field is constant over 
the detection length and taking the quantity of (f /Io)/V\\ to be 
the mean value of the measurements for oscillations 1 and 2, we 
obtain = 59° + 8° from Eq. ©. Then with Eq. © the sound 
speed is estimated to be 128+25 km s"^ From Cs = 0.1527^^^ 
km s"^ we estimate the average temperature of the loop near the 
footpoint about 0.7+0.3 MK. This result is in good agreement 
with what is obtained from observa tions of 3D propag ation of 
similar waves with STEREO/EUVI (iMarsh et al. 1120091) . 

4. Discussion 

We have reported on the first observation of the multiple-period 
(12 and 25 min) propagating disturbances along a fan-like coro- 
nal structure in both intensity and Doppler shift in the Fe xii 195 
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Fig. 3. Wavelet analysis for oscillation 1 at 3;= 150 '. (a) Time 
profiles of Doppler shift and relative intensity. Here positive val- 
ues for the Doppler shift represent the blueshifted emission, (b) 
The wavelet power spectrum for Doppler shift. The dark color 
represents high power and the dotted contour encloses regions 
of greater than 99% confidence for a white-noise process, (c) 
The global wavelet spectrum (solid line) and its 99% confidence 
level (dotted line), (d) and (e): The same as (b) and (c) but for 
relative intensity. 



A line with Hinode/EIS. The phase and amplitude relationships 
between the intensity and Doppler shift oscillations provide con- 
vincing evidence that they are a propagating slow magnetoa- 
coustic wave. 

With high cade n ce (15 s) EIT 195 A observations, 
iBerghmans & Clettel (1199 9) first detected the propagating in- 
tensity disturbances in fanlike coronal loops. They noticed that 
these disturbances are associated with recurrent footpoint bright- 
enings, appearing as transient flows ejected periodically. In our 
case the propagating features are detected in both intensity and 
Doppler shift. The small Doppler shift amplitude of several 
km s"^ is consistent with the intensity amplitude as expected for 
a slow magnetoacoustic wave, thus, interpreting these oscilla- 
tions in terms of propagating slow magnetosonic waves is more 
appropriate. In addition, the apparent propagation speed of the 
oscillations >100 km s"^ is an order of magnitude higher than 
any flow speeds detected in this region. 

In our case the dominant oscillations have periods of about 
25 min, much lon ger than those observed in the TRACE loops 
of around 5 min (IDe Moortel et al. ll2Q02h . It is not clear why 
the longer periods rep orted here were no t observed in TRACE 
loops. One reason for iDe Moortel et al. I (|2b02) detecting only 
the 5 min oscillations could be attributed to the selection crite- 
rion they used, where a running diff'erence technique with time 
lag of 90 s would elimi nate the long-period co mponents that 
may exist in some cases. De Pontieu et al. I (l2005l) suggest that 5 
min waves in coronal loops come from the leakage of /?-modes 
through the chromosphere tunneled by highly inclined magnetic 
fields. However, the origin of 20-30 min waves is difficult to ex- 
plain by such a mechanism because their period is much longer 
than the cutoff' period of the slow magnetoacoustic waves in the 
chromosphere (Bel & Leroy 1977). Therefore, another explana- 
tion is required, e.g., possibly the coupling with the photospheric 
/-mode caused by local convective motions. 
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Fig. 4. Wavelet analysis for oscillation 2 at y=110'\ The anno- 
tations are same as in Fig.O 



The detection length is deffned as the distance along the 
loop over which the intensity disturbances can be observed. 
TRACE observations show that the waves are damped quickly 
with a detection length typically in the range of 3-23 Mm 
( iDe Moortel ^TZI [2QQ2h . Theoretical and numerical studies 
show that thermal conduction can play an important role in 
the quick damping of slo^v magnetosonic waves in coronal 
loops (q.2. [Ofman & Wanj] I2OO2I: iDe Moortel & HoodI 120041: 
iKlimchuk et al. |[2004|) . In our case the detection length is mea- 
sured to be about 69 Mm for intensity and about 87 Mm for 
Doppler velocity, much longer than in the TRACE case. This 
diff'erence may be explained by the dependence of the damp- 
ing lengt h on the wavelength an d the temperature of the ambient 
plasma ( Klimchuk et al. 1120041) . Thus, the longer period waves 
reported here have longer wavelength than in previous obser- 
vations, resulting in weaker damping by thermal conduction. In 
addition, the present loop appears cooler (~ 0.7 MK as estimated 
in the previous section) than previously observed loops (1.0-1.5 
MK for TRACE loops). 
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Fig. A.l. Comparisons between EIS and TRACE observations, 
(a) Constructed TRACE 195 A time-distance image along the 
EIS slit, (b) Same as (a) but with the background subtraction, 
where the background is taken as the average over time at each 
position along the cut. (c) The EIS Fe xii intensity time series, 
(d) Same as (c) but with the background subtraction. 
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Fig.A.2. Wavelet analysis for oscillation 3 averaged over 
y=60 "-65 ". The annotations are same as in Fig.O 
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a) Doppler velocity and intensity time series 
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Fig. A.3. Wavelet analysis for oscillation 1 at y=150'', which is 
the same as Fig.[3]but with a subtraction of the background taken 
as a 30 min average smoothing. 

Table A.l. Measurements of the oscillations in fanlike loops 
found in the EIS sit-and-stare data in the Fe xii (195 A) line^. 



Oscillation 


Pv 


Pi 


Ay 


A/ 


No. 


(min) 


(min) 


(km s-i) 


(%) 


1 


25.6+6.3 


26.2+4.9 


1.6+0.5 


4.9+1.5 


2^ 


25.1+7.0 


26.2+5.0 


1.3+0.3 


3.8+0.6 




11.9+2.8 


13.5+4.7 






3 


11.2+3.8 


11.2+2.9 


1.0+0.3 


2.8+0.7 



^ Py and Pi are periods and Ay and A/ are amplitudes for the 

Doppler shift and relative intensity oscillations, respectively. 
^ Two harmonic components are found for oscillation 2. 



